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The present work shows a novel fractal dimension method for shape analysis. The 
proposed technique extracts descriptors from the shape by applying a multiscale ap- 
proach to the calculus of the fractal dimension of that shape. The fractal dimension 
is obtained by the application of the curvature scale-space technique to the original 
shape. Through the application of a multiscale transform to the dimension calculus, 
it is obtained a set of numbers (descriptors) capable of describing with a high preci- 
sion the shape in analysis. The obtained descriptors are validated in a classification 
process. The results demonstrate that the novel technique provides descriptors highly 
reliable, confirming the precision of the proposed method. 
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I. INTRODUCTION 



One of the most important research areas in pattern recognition and image analysis is 
the analysis of shapes. Shapes are capable of describing an object by preserving the most 
relevant information about it, simplifying its aspect. This enables us to use computational 
algorithms to solve important problems in computational visiorP. 

Particularly, the literature presents various methods aiming to provide a single number 
or a set of them, which can be used to describe the essential characteristics of the shape. 
These numbers are called shape descriptors. Among the methods for the obtainment of 
these descriptors, one that has a growing number of works showed in the literature is the 
fractal descriptor^. 

Formally defined irP, the fractal geometry studies mathematical entities, which cannot 
be perfectly modelled and described by the conventional Euclidean geometry. In the Nature, 
we can find a number of these objects, such as the sky clouds, the trees, the river streams, 
etc. MandelbrotP presents fractal objects as complex structures characterized mainly by the 
self-similarity, i.e., each part of the object is a copy in reduced scale of the original object. 
In authentic fractal objects, the self-similarity property is seen in infinite levels of scale. 

In his study about fractal geometry, MandelbrotP showed more formally the concept of 
fractal dimension and the possibility of its use as a measure to describe fractal objects. IrP, 
it is still demonstrated that the fractal geometry and, consequently, the fractal dimension 
concept, could be extended to objects that are not necessarily true fractal entities, but, 
which demonstrate characteristics of fractals as the self-similarity (although to a limited 
degree). With the advent of computers and the representation of natural objects by digital 
images, the study of fractal methods provided an important source of descriptors for objects 
in natural scenes. Recently, we can find several works in the literature using the fractal 
dimension as a descriptor of objects in applications at different areas of the scienc^^. 

In the specific case of shape analysis,^ presents a novel approach of fractal descriptors, 
where the concept of Multi-scale Fractal Dimension (MFD) is defined. In that approach, 
the fractal dimension is not used as a unique descriptor of an object but the dimension is 
calculated at different levels of observation (from a micro to a macroscopic scale). Thus, 
a set of numbers (each one corresponding to a different observation scale) is generated to 
compose a vector of fractal descriptors. This approach has demonstrated to represent the 
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shape with more richness than the use of a simple fractal dimension descriptor and many 
other works have followed the same approach and showed interesting results in pattern 
recognition problem j^ * ^ * ^ l 

The present work proposes a novel MFD approach based on the fractal dimension calculus 
method developed inP. InP the fractal dimension is calculated by the multi-scale approach 
called Curvature Space Scale (CSS)I3S3, in which the fractal dimension is calculated by the 
determination of the value of geometric curvature of the shape at a varying scale of obser- 
vation. The proposal of this work is to use the set of curvature measures as a set of MFD 
descriptors for the shape. The method is validated in a task of classification of object shapes 
over a dataset whose classes are previously known. 

The next sections show the main mathematical concepts used in this work, i.e., curvature 
and fractal dimension. It also describes the CSS technique and how it can be used for 
the calculus of fractal dimension. The following section describes the Multi-scale Fractal 
Dimension method. An experiment to validate our approach is set on the posterior section. 
The next section presents the results and discussion and the last one does the conclusions 
of the work. 

II. SHAPE AND FRACTAL DIMENSION 

The use of shapes to characterize objects in a scene constitutes an important branch of 
pattern recognitiorP. The shape of an object can be described as a binary image representing 
the silhouette of the analyzed object. It contains the main information relative to the original 
object and shows less complexity than the representation in the original image. 

In pattern recognition applications, it is interesting to represent a shape by a set of nu- 
merical descriptors, i.e., a feature vector. The literature shows various techniques to generate 
these descriptors. Particularly, one interesting shape descriptor is the fractal dimension. 

The concept of fractal objects dates back to the 19*'^ century. However, MandelbroiP is 
the first known author to show a more formally detailed study. Fractals are objects with 
infinite complexity, generated from simple Euclidean geometric entities and following, in 
most cases, a simple rule of construction. They are also characterized by their self-similarity, 
i.e., if one takes a small piece of the object, it can notice that this piece is an approximated 
copy, reduced in scale, of the original object. In true fractal objects, this property can be 
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observed at all levels of scale (infinitely). 

Mandelbrot noticed in his work that the fractals were an interesting tool to model objects 
found in the nature. In fact, we can easily observe that there are various objects around 
us, which are composed by the simple repetition of a pattern. Thus, these objects are, 
in a certain sense, similar to fractal objects, presenting the self-similarity property. It is 
important to notice that the natural objects cannot be considered true fractals once its 
self-similarity is manifested only at a limited number of scales. They are pseudo-fractal 
objects.l^ demonstrates the efficiency of the modelling of objects of the real world by tools 
of the fractal geometry. 

In problems of computational vision, an important task is the identification and classifi- 
cation of objects in a natural scene. In the conventional process of a computational vision 
system, an important step is the calculus of numbers capable of describing the analyzed ob- 
ject, the descriptors. When the object is modelled as a pseudo-fractal entity, an important 
descriptor that can be immediately obtained is the fractal dimension. The concept of fractal 
dimension was defined for true fractal objects and is formally defined irP as: 

D = lim (1) 

where L is the length of the fractal object, u is the length of a measure line smaller than 
L and is the number of objects of length u necessary to "cover" the fractal object. In the 
case of pseudo-fractal objects the literature shows the use of some methods based on the 
original concept, as the Box-countin^^^, Minkowski mass-radiu^, and others. 

In practice, fractal dimension demonstrates to be a measure of the spatial occupation 
of the object or still of its complexity. Any method used to estimate the fractal dimension 
must capture the complexity of the shape at different levels of observation (micro and macro 
scales). Such process elucidates the importance of the fractal dimension in the description 
of the geometric structure of the shape, consequently, distinguishing it from other objects, 
as required by the pattern recognition algorithms. 



A. Curvature Scale Space 



Curvature is a powerful tool in the shape analysis field. Its use enables us to study the 
behavior of a curve by its changes in orientation. Thus, it is possible to detect relevant 
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points in the shape which can be used for its description and identification. The curvature 
is defined in terms of derivatives as 

(x«(t)2+y(l)(t)2)3/2 ' 

where "-^^ and denote the first and second derivatives, respectively, and C(t) = {x(t),y(t)) 
is a parametric vector representing the curve under analysis. 

The study of the curvature property in the scale space has shown great robustness in the 
identification of shapes under many types of transformations, such as scaling, orientation 
changes, translation and even noise. The Curvature Scale Space (CSS) is achieved by apply- 
ing a multi-scale transformation over the original curve. This transformation is, basically, 
the convolution of a Gaussian function g{t, a) over the parametric curve. This is performed 
in order to reduce the effects of noise and high frequency information before the curvature 
measurementP^*^^*^ 



X{t,a) = x{t)*g{t,a), (3) 
Y{t,a)=yit)*g{t,a), (4) 

where 

1 

^(^' ^) = — As= exp^ , (5) 

(JV/TT 

is the Gaussian function with standard deviation a and * denotes convolution. 
Thus, curvature equation considering the scale evolution is defined as 

, XW(t,a)y(^)(t,a)-X(^)(t,a)yW(t,a) 
^ ' ^ (X(i)(t,a)2 + y(i)(t,a)2)3/2 ' W 

where the value of a is linearly increased to achieve the set of curves, which composes 

the Curvature Scale Space of a given shape (Figure [T|. 



B. Estimating fractal dimension by Curvature 

In the CSS approach, a Gaussian function is applied over the parametric curve in order to 
reduce the amount of high frequency information and noise present in the original shape. As 
the value of the a parameter in the Gaussian function decreases, the higher is the smoothing 
level applied over the curve. During this process, the aspect of the shape changes. We note 
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FIG. 1. (a) Fish coiitours;(c) Curvature Scale Space. 

that the higher changes in the contour orientation are diminished due to the reduction of 
the information contained in the curve. As a result, the shape becomes more similar to a 
circle (Figure [2]). This process reflects in the complexity of the shape, so that, curvature 
scale space can be used to estimate the fractal dimension of a curve. 




(a) (b) (c) (d) 

FIG. 2. Smoothing of a shape contour as the a parameter changes. 

In the Bouligand-Minkowski method, the shape of an object is dilated and its area com- 
puted for different dilatation radius r to estimate the shape complexity (Figure |3]). A similar 
idea can be applied to the curvature. Like in the Bouligand-Minkowski method, we define 
the amount of curvature S{a) as the sum of the module of all curvature coefficients computed 
for a shape: 
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S{a) = Y.\k{t,cr)\, (7) 

t 

where S{a) represents the resuhing orientation of the curve at scale a. For each a value 
used, a different amount of curvature S{a) is computed. However, defining a set of a values 
to be used is a difficult task. Thus, as in CSS method, we opt to linearly increase its value. 

By accumulating the subsequent values for S'(cr), it is possible to create a relation of 
dependence among scales, where posterior scales are influenced by previous scales. This 
is just similar to the process of dilating a shape performed by the Bouligand-Minkowski 
method, where a dilatation using a radius r is the accumulation of every point at a distance 
from the shape ranging from to r. Thus, the accumulated amount of curvature is defined 
as 



Aia) = Y.Si^). (8) 

i=0 

The degree of space filling of a curve can be quantified by considering its order of growth. 
For this, it is necessary to analyze the behavior of o" x A{a) curve in the logarithmic scale. 
Figure |4] shows the log-log plot of the a x A(a) curve computed for the Koch snowflake (or 
Koch star). The Koch snowflake is a mathematical curve of infinite length and one of the 
earliest fractal curves described in the literature. At each iteration, one line segment of this 
curve is replaced by four line segments, each one having one-third of the size of its previous 
stage, what leads to the fractal dimension d = log4/log3 = 1.2618. Note that the curve 
achieved obeys a power law. According to^^, this power law can be used to estimate the 
fractal dimension of the original shape: 

d = hm (9) 
f^^o logo- ^ 

By applying a line regression over the log-log plot of cr x A{a), it is possible to estimate 
a straight line with slope a, where d = a is the estimated fractal dimension for the closed 
curve C{t). Notice that the method estimated the fractal dimension with d = 1.2644. This 
result is very close to the theoretical fractal dimension of the Koch snowflake and shows the 
accuracy of the proposed method and its straight relation to the fractal theory. 
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(c) 



(d) 



FIG. 3. Process of shape dilation by the BouHgand-Minkowski method (Hght gray, points corre- 
spondent to the current radius r; dark gray, accumulated area): (a) Original shape (r = 0); (b) 
r2 = 1; (b) = 2; (b) = 4. 




(a) (b) 
FIG. 4. (a) Koch snowflake [d = 1.2618); (b) Curvature Fractal Dimension [d = 1.2644). 

III. MULTI-SCALE FRACTAL DIMENSION 



As previously described, the changes in the value of the a parameter in the Gaussian 
function affects the aspect of the shape (Figure |2]). As a result, small variations are incor- 
porated to the log-log curve a x A{cr). Thus, eventually, a single non-integer number d may 
be not adequate to characterize all the complexity richness existent in an object. Moreover, 
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non-fractal objects, such as a shape curve, has finite size. This imphes that its dimension 
goes to zero as the visuahzation scales increase. This is corroborated by the tendency of 
the shape to become more similar to a circle as the the information of high frequencies is 
reduced. 

In order to solve this deficiency of the fractal dimension methods and aiming to provide 
a better description of objects in terms of complexity, the Multi-Scale Fractal Dimension 
has been proposecP^^'^. This approach exploits the infinitesimal limit in the linear interpo- 
lation by using the derivative, thus yielding a function that binds the changes in the object 
complexity to the visualization scale changes (Figure |5]). This function enables us to perform 
a more effective discrimination of the object and it is defined as: 



d\ogA{a) 



dloga 

where d{a) represents the complexity of the object at scale a. 



0,5 

loser 



1 1.5 2 2.5 




log a 



(a) (b) 
FIG. 5. (a) log-log curve, (b) Multi-scale Fractal Dimension. 



IV. EXPERIMENTS 

An experiment was performed to evaluate the proposed fractal approach as a feasible 
shape descriptor. The experiment aimed the classification of a large database composed 
by shapes under rotation and scale transformations. The database consists of 1100 classes, 
where each class is one of the 1100 fish contours obtained from the database available atP^. 
Each class is composed of 10 different manifestations (5 rotations and 5 scales) of the fish 
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contour associated, thus resulting in a database containing 11000 shapes. Figure [6] shows 
some examples of fish contours present in the database. 

Statistical evaluation of the Multi-scale curves were carried out using the Linear Discrim- 
inant Analysis (LDA). The LDA is a supervised statistical classification methocP^'^, which 
searches a linear combination of the descriptors (independent variables) that results in its 
class (dependent variable). Thus, the LDA method attributes an observation x to the class 
i which presents the higher conditional probability: 

p{t\x)>p{j\x),yj^t. (10) 

where P{i\x) is the conditional probability of i given x. Considering that P{i\x) obeys 
a multivariate gaussian distribution with a single covariance matrix for all the classes, an 
object is attributed to the class i which provides the higher value for the function ff 

fi = pi,C~^xl - O.S/iiC-^Xfc + log(pi), (11) 

where /ij is the mean of the descriptors of class i, C is the covariance matrix of the data 
set and pi is the a priori probability of class i. 

The statistical analysis was carried out over the samples using a K-fold cross-validation 
scheme with 5 folds. In order to provide a fair evaluation of the properties of the proposed 
method (such as the tolerance to rotation and scaling), one manifestation of rotation and 
scaling was considered in each class, in each fold. 




FIG. 6. Examples of fish images used in the experiments. 
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V. RESULTS AND DISCUSSION 



To provide a robust evaluation of the proposed approach, a comparison with the 
Bouhgand-Minkowski method was considere d^ * ^^ * ^^ l The Bouhgand-Minkowski is one of 
the most accurate methods to estimate the fractal dimension. In this method, the shape 
of an object is dilated by a disc of radius r and the influence area is used to estimate the 
shape complexity. The influence area is very sensitive to structural changes of the shape, 
thus even small changes in the shape structure can be detected. For this experiment, we 
considered a dilation radius r = 50. 

Finite Difference method-^ was applied over log-log curves computed for both Bouhgand- 
Minkowski and the proposed approach, thus resulting in a Multi-scale fractal dimension 
curve characteristic for that shape. A total of 50 descriptors was considered from Multi- 
scale curves for shape characterization. 

Table |T] presents the classification results yielded for each approach. Variations in the 
size or orientation of a shape are a very common issue in shape analysis. Tolerance to these 
transformations is a fundamental characteristic for a method to be considered robust in 
shape analysis applications. As we can see, the Bouhgand-Minkowski method is not scale 
invariant. Thus, in order to reduce this problem and to perform a fairer comparison, a 
second experiment was performed where the shapes were first normalized according to their 
diametei^. 

TABLE I. Comparison results for different complexity estimation methods. 

Method Samples correctly classified Success rate 

Proposed approach 10685 97.14 

Bouligand-Minkowski 1573 14.30 

Bouligand- Minkowski (normalized) 8807 80.06 

Both methods compared are tolerant to orientation changes in the shape. Nevertheless, 
it is important to emphasize that images exist only in the discrete space, and they cannot 
be freely rotated. A small error is always added to the multi-scale curve due to shape 
rotation. The tolerance in the Bouligand-Minkowski method is due to the use the Euclidean 
distance during the computing of the influence area. This is a distance measure that is not 
affected by transformations in orientation. In our proposed approach, the rotation tolerance 
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is most since curvature is not sensitive to orientation changes in the shape. The curvature 
is a measure related to the curve orientation at a specific point, i.e., the changes in the 
direction inside the curve. Thus, the changes in the orientation of the shape do not affect 
the orientation among its inner segments. 

Unhke the Bouhgand-Minkowski method, which needs a normahzation step performed 
directly over the shape image, the proposed approach is tolerant to scale changes and so no 
normalization of the data is necessary. The main reason for its scaling tolerance is the use 
of the sum of the module of the curvature coefficients as the basis to estimate the shape 
complexity (Figure [?]). As in the rotation tolerance, scaling changes does not affect the 
orientations of the inner structures of the contour. It only increases when the size of the 
shape increases, i.e., the number of points in its contour. Minimum and maximum local 
points are also preserved when the direction changes in the contour as well as the sum of its 
coefficients. Thus, the behavior of the multi-scale fractal dimension curve is preserved even 
for scaled shapes. 




FIG. 7. Curvature and the sum of its coefficients, S{a), for two contour at different scales. 
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VI. CONCLUSION 



This work presented the study of a novel technique for extracting shape descriptors using 
a fractal methods based on the curvature scale space concept. 

Results demonstrated that the proposed technique achieved a great accuracy in the classi- 
fication of the shape database containing rotated and scaled samples, showing an interesting 
correctness rate of 97.14%. The method also showed the best accuracy in comparison to 
the Bouligand-Minkowski method, as it is being intrinsically invariant to scale, a property 
which cannot be found in the conventional Bouligand-Minkowski method. 

As the Bouligand-Minkowski method is shown in the literature as an efficient technique to 
obtain fractal shape descriptors, the results presented in this work suggest that the proposed 
technique is a worthy option for providing shape descriptors for classification tasks, as it is 
invariant to rotation and scale, operations commonly found in practical situations. Such fact 
turns possible the application of the proposed technique to various problems in computer 
vision involving the classification of shapes. 

We still suggest, as a future work, a deeper study about this novel method, in order to 
verify its application to another pattern recognition and computer vision problems. 

ACKNOWLEDGMENTS 

A.R.B. acknowledges support from FAPESP (2006/54367-9). J.B.F. acknowledges sup- 
port from CNPq (870336/1997-5) and FAPESP (2006 / 53959-0). O.M.B. acknowledges 
support from CNPq (Grant #308449/2010-0 and #473893/2010-0) and FAPESP (Grant # 
2011/01523-1). 

REFERENCES 

^S. Loncaric, "A survey of shape analysis techniques," Pattern Recognition, 31, 983-1001 
(1998). 

2R. de O. Plotze, M. Falvo, J. G. Pdua, L. C. Bernacci, M. L. C. Vieira, G. C. X. Ohveira, 
and O. M. Bruno, "Leaf shape analysis using the multiscale minkowski fractal dimension, 
a new morphometric method: a study with passifiora (passifioraceae)," Canadian Journal 
of Botany, 83, 287-301 (2005). 

13 



^B. B. Mandelbrot, The Fractal Geometry of Nature (Freeman, New York, 1968) p. 468. 
^M. Carlin, "Measuring the complexity of non-fractal shapes by a fractal method," PRL: 
Pattern Recognition Letters, 21, 1013-1017 (2000). 

^L. Billeci, G. Pioggia, F. Vaglini, and A. Ahluwalia, "Assessment and comparison of neural 
morphology through metrical feature extraction and analysis in neuron and neuron-glia 
cultures." Journal of biological physics, 35, 447-464 (2009). 

6J.-E. Qu, X.-P. Guo, Z.-Y. Chen, G.-W. Song, and X.-H. Zhang, "Fractal research on 
afm images of polycrystalline aluminum surface with adsorption film of sds inhibitor," 
Materials Science, 45, 114-124 (2009). 

'''R. da S. Torres, A. F. ao, and L. da F. Costa, "A graph-based approach for multiscale 
shape analysis," Pattern Recognition, 37, 1163-1174 (2003). 

^A. R. Backes and O. M. Bruno, "Medical image retrieval based on complexity analysis (to 
appear)," Machine Vision and Applications, 21, 217-227 (2010). 

^A. R. Backes, J. B. Florindo, and O. M. Bruno, "A novel approach to estimate fractal 

dimension from closed curves," in CAIP (2009) pp. 253-260. 
^°L. da F. Costa and R. M. C. Jr., Shape Analysis and Classification: Theory and Practice 

(CRC Press, Boca Raton-FL, 2000). 
^^C. Tricot, Curves and Fractal Dimension (Springer- Verlag, New York, 1995). 
^^A. P. Witkin, "Scale space filtering: a new approach to multi-scale descriptions," in 

ICASSP - IEEE International Conference on Acoustics, Speech, and Signal Processing 

(GRETSI, 2003) pp. 79-95. 
^^F. Mokhtarian and S. Abbasi, "Matching shapes with self-intersections: application to leaf 

classification," IEEE Transactions on Image Processing, 13, 653-661 (2004). 
^^C. W. Emerson, N. N. Lam, and D. A. Quattrochi, "Multi-scale fractal analysis of im- 
age texture and patterns," Photogrammetric Engineering and Remote Sensing, 65, 51-62 

(1999). 

^^R. C. Gonzalez and R. E. Woods, Digital Image Processing, 2nd ed. (Prentic-Hall, New 
Jersey, 2002). 

"http: / / www.ee.surrey.ac.uk/cvssp / demos/css/demo.html," (2003) . 
-"^^B. S. Everitt and G. Dunn, Applied Multivariate Analysis, 2nd ed. (Arnold, London, 2001). 
^^K. Fukunaga, Introduction to Statistical Pattern Recognition, 2nd ed. (Academic Press, 

San Diego, 1990). 



14 



'O. M. Bruno, R. de O. Plotze, M. Falvo, and M. de Castro, "Fractal dimension applied 
to plant identification," Information Sciences, 178, 2722-2733 (2008). 
'G. D. Smith, Numerical Solution of Partial Differential Equations: Finite Difference Meth- 
ods, 3rd ed. (Oxford Press, Oxford-USA, 1986). 



15 



